Understanding Original Antigenic Sin in Influenza with a Dynamical System

Original antigenic sin is the phenomenon in which prior exposure to an antigen leads to a subsequent suboptimal immune response to a related antigen. Immune memory normally allows for an improved and rapid response to antigens previously seen and is the mechanism by which vaccination works. I here develop a dynamical system model of the mechanism of original antigenic sin in influenza, clarifying and explaining the detailed spin-glass treatment of original antigenic sin. The dynamical system describes the viral load, the quantities of healthy and infected epithelial cells, the concentrations of naïve and memory antibodies, and the affinities of naïve and memory antibodies. I give explicit correspondences between the microscopic variables of the spin-glass model and those of the present dynamical system model. The dynamical system model reproduces the phenomenon of original antigenic sin and describes how a competition between different types of B cells compromises the overall effect of immune response. I illustrate the competition between the naïve and the memory antibodies as a function of the antigenic distance between the initial and subsequent antigens. The suboptimal immune response caused by original antigenic sin is observed when the host is exposed to an antigen which has intermediate antigenic distance to a second antigen previously recognized by the host's immune system

Computational and Theoretical Analysis of Influenza Virus Evolution and Immune System Dynamics

Influenza causes annual global epidemics and severe morbidity and mortality. The influenza virus evolves to escape from immune system antibodies that bind to it. The immune system produces influenza virus specific antibodies by VDJ recombination and somatic hypermutation. In this dissertation, we analyze the mechanism of influenza virus evolution and immune system dynamics using theoretical modeling and computational simulation. The first half of this thesis discusses influenza virus evolution. The epidemiological data inspires a novel sequence-based antigenic distance measure for subtypes H1N1 and H3N2 virus, which are superior to the conventional measure using hemagglutination inhibition assay. Historical influenza sequences show that the selective pressure increases charge in immunodominant epitopes of the H3 hemagglutinin influenza protein. Statistical mechanics and high-performance computing technology predict fixation tendencies of the H3N2 influenza virus by free energy calculation. We introduce the notion of entropy from physics and informatics to identify the epitope regions of H1-subtype influenza A with application to vaccine efficacy. We also use entropy to quantify selection and diversity in viruses with application to the hemagglutinin of H3N2 influenza. Using the bacterial E. coli as a model, we show the evidence for recombination contributing to the evolution of extended spectrum β-lactamases (ES-BLs) in clinical isolates. A guinea pig experiment supports the discussion on influenza virus evolution. The second half of the thesis discusses immune system dynamics. We design a two-scale model to describe correlation in B cell VDJ usage of zebrafish. We also introduce a dynamical system to model original antigenic sin in influenza. This dissertation aims to help researchers understand the interaction between influenza virus and the immune system with a quantitative approach

Demand Driven Research Support in the Library: A Case Study in History and Digital Humanities

This section presents the timeline of instruction support, research consultation, and collection development offered by faculty librarians Molly Castro and Christopher M. Jimenez to support History Professor Kyle Pan and History grad student Noel Hernandez. Noel’s dissertation research and accompanying digital history project on the history of rock bands in Latin America serve as a case study on demand driven research support from librarians in both Information & Research Services and Digital Humanities

A Multi-Scale Model for Correlation in B Cell VDJ Usage of Zebrafish

The zebrafish (\emph{Danio rerio}) is one of the model animals for study of immunology because the dynamics in the adaptive immune system of zebrafish are similar to that in higher animals. In this work, we built a multi-scale model to simulate the dynamics of B cells in the primary and secondary immune responses of zebrafish. We use this model to explain the reported correlation between VDJ usage of B cell repertoires in individual zebrafish. We use a delay ordinary differential equation (ODE) system to model the immune responses in the 6-month lifespan of a zebrafish. This mean field theory gives the number of high affinity B cells as a function of time during an infection. The sequences of those B cells are then taken from a distribution calculated by a "microscopic" random energy model. This generalized $NK$ model shows that mature B cells specific to one antigen largely possess a single VDJ recombination. The model allows first-principles calculation of the probability, $p$, that two zebrafish responding to the same antigen will select the same VDJ recombination. This probability $p$ increases with the B cell population size and the B cell selection intensity. The probability $p$ decreases with the B cell hypermutation rate. The multi-scale model predicts correlations in the immune system of the zebrafish that are highly similar to that from experiment.Comment: 29 pages, 10 figures, 1 tabl

Lattice-free prediction of three-dimensional structure of programmed DNA assemblies

DNA can be programmed to self-assemble into high molecular weight 3D assemblies with precise nanometer-scale structural features. Although numerous sequence design strategies exist to realize these assemblies in solution, there is currently no computational framework to predict their 3D structures on the basis of programmed underlying multi-way junction topologies constrained by DNA duplexes. Here, we introduce such an approach and apply it to assemblies designed using the canonical immobile four-way junction. The procedure is used to predict the 3D structure of high molecular weight planar and spherical ring-like origami objects, a tile-based sheet-like ribbon, and a 3D crystalline tensegrity motif, in quantitative agreement with experiments. Our framework provides a new approach to predict programmed nucleic acid 3D structure on the basis of prescribed secondary structure motifs, with possible application to the design of such assemblies for use in biomolecular and materials science.United States. Office of Naval Research (ONR N000141210621)National Science Foundation (U.S.) (NSF-DMREF Program CMMI1334109

Selective Pressure to Increase Charge in Immunodominant Epitopes of the H3 Hemagglutinin Influenza Protein

The evolutionary speed and the consequent immune escape of H3N2 influenza A virus make it an interesting evolutionary system. Charged amino acid residues are often significant contributors to the free energy of binding for protein–protein interactions, including antibody–antigen binding and ligand–receptor binding. We used Markov chain theory and maximum likelihood estimation to model the evolution of the number of charged amino acids on the dominant epitope in the hemagglutinin protein of circulating H3N2 virus strains. The number of charged amino acids increased in the dominant epitope B of the H3N2 virus since introduction in humans in 1968. When epitope A became dominant in 1989, the number of charged amino acids increased in epitope A and decreased in epitope B. Interestingly, the number of charged residues in the dominant epitope of the dominant circulating strain is never fewer than that in the vaccine strain. We propose these results indicate selective pressure for charged amino acids that increase the affinity of the virus epitope for water and decrease the affinity for host antibodies. The standard PAM model of generic protein evolution is unable to capture these trends. The reduced alphabet Markov model (RAMM) model we introduce captures the increased selective pressure for charged amino acids in the dominant epitope of hemagglutinin of H3N2 influenza (R2 > 0.98 between 1968 and 1988). The RAMM model calibrated to historical H3N2 influenza virus evolution in humans fit well to the H3N2/Wyoming virus evolution data from Guinea pig animal model studies

Evolution of H3N2 Influenza Virus in a Guinea Pig Model

Studies of influenza virus evolution under controlled experimental conditions can provide a better understanding of the consequences of evolutionary processes with and without immunological pressure. Characterization of evolved strains assists in the development of predictive algorithms for both the selection of subtypes represented in the seasonal influenza vaccine and the design of novel immune refocused vaccines. To obtain data on the evolution of influenza in a controlled setting, naïve and immunized Guinea pigs were infected with influenza A/Wyoming/2003 (H3N2). Virus progeny from nasal wash samples were assessed for variation in the dominant and other epitopes by sequencing the hemagglutinin (HA) gene to quantify evolutionary changes. Viral RNA from the nasal washes from infection of naïve and immune animals contained 6% and 24.5% HA variant sequences, respectively. Analysis of mutations relative to antigenic epitopes indicated that adaptive immunity played a key role in virus evolution. HA mutations in immunized animals were associated with loss of glycosylation and changes in charge and hydrophobicity in and near residues within known epitopes. Four regions of HA-1 (75–85, 125–135, 165–170, 225–230) contained residues of highest variability. These sites are adjacent to or within known epitopes and appear to play an important role in antigenic variation. Recognition of the role of these sites during evolution will lead to a better understanding of the nature of evolution which help in the prediction of future strains for selection of seasonal vaccines and the design of novel vaccines intended to stimulated broadened cross-reactive protection to conserved sites outside of dominant epitopes